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We study the dominant effect of a long wavelength density perturbation S(Xl) on short distance 
physics. In the non-relativistic. limit, the result is a uniform acceleration, fixed by the equivalence 
principle, and typically has no effect on statistical averages due to translational invariance. This same 
reasoning has been formalized to obtain a “consistency condition” on the cosmological correlation 
functions. In the presence of a feature, such as the acoustic peak at Ibao, this naive expectation 
breaks down for A l < £bao- We calculate a universal piece of the three-point correlation function in 
this regime. The same effect is shown to underlie the spread of the acoustic peak, and is calculable to 
all orders in the long modes. This can be used to improve the result of perturbative calculations — a 
technique known as “infra-red resummation”— and is explicitly applied to the one-loop calculation 
of power spectrum. Finally, the success of BAO reconstruction schemes is argued to be another 
empirical evidence for the validity of the results. 


I. INTRODUCTION 

Local experiments performed in a small laboratory 
cannot reveal the existence of the uniform gravitational 
field of a long-wavelength matter density perturbation, 
e.g. Sl( x, t) = 8 q (t)cos(q ■ x). By the equivalence prin¬ 
ciple, the laboratory and all of its belongings fall with 
a uniform acceleration — V$ i (a:i ab , t), where ^ B {x,t) = 
—^ttGci 2 p(t)SL(x, t)/q 2 , and p(t) is the mean matter den¬ 
sity of the Universe. However, two distant laboratories 
with separation larger than 1/q experience different ac¬ 
celerations. A distant observer sees a clear correlation 
between the relative motion of the two and the underly¬ 
ing density perturbation. 

The motion in the field of a long-wavelength mode is 
easiest to find from the fact that everything falls in the 
same way as a dark matter particle does. Possible devia¬ 
tions are suppressed by additional derivatives of the long 
mode. For dark matter, the linearized continuity equa¬ 
tion implies v ~ ~ y&A The total displacement since 
t = 0 is then 

Ax = 5 q (t) sin(g • x) q/q 2 . (1) 

The small laboratories of the cosmologist, like stars 
and galaxies, are observed at a single point in their lifes¬ 
pan. Hence, the relative motion of any given pair is im¬ 
possible to determine. What is possible is to see how the 
distribution of pairs is correlated with Sl • For pairs of 
any objects, say galaxies, equation (1) implies 

-&(*>*) 

( 2 ) 

+2(5, (t) sin ^ • V£ s (x,t), 

where £ g (x,t) is an average 2-point correlation function. 
Not surprisingly, the distribution of pairs with separa¬ 
tion much less than the long wavelength, q ■ x <C 1, is 
hardly effected by the long mode. The second line would 
in this case correspond to the effect of living in an over 
(under) dense Universe. An effect of order <$z,x|V£ g |, 


which for an approximately scale invariant spectrum, 
|V£ g (£e, t) | ~ £ g (x, t)/x, is comparable to dynamical con¬ 
tributions of order d^g, which are neglected anyway on 
the right-hand side (r.h.s.). However, even if q ■ x 1, 
when we do expect the long-wavelength mode to induce 
a large relative motion, the second line of (2) is often 
negligibly small. Scale invariance now implies that it is 
of order S^g/qx. 

The relative motion is noticeable only if the distribu¬ 
tion of pairs has a feature such that the derivative in the 
second line of (2) becomes large. One such feature does 
exist in the Universe at the baryon acoustic oscillation 
(BAO) peak. For x ~ ^bao, 

|V£ 9 | ~ Ug » (3) 

o £bao 


where a is the width of the peak. At this separation, 
the effect of the long mode on the distribution of pairs 
is of order 5 L £ B AO^g/^ for q ^baO’ anc ^ ^L^glo® fo r 
^bao Q •C cr — 1 , which are both dominant compared 
to the 0(6L^,g) dynamical effects. In what follows, we 
explore the implications of this simple observation for 
the shape of the correlation functions around the BAO 
scale, and its connection to broadening of the peak. 1 


II. CORRELATION WITH THE LONG MODE 

Real space. — An approximate three-point correlation 
function can be obtained in this regime by correlating (2) 


1 The initial time for this problem can be taken long after the 
recombination, when the acoustic peak is already in place, but 
the modes of interest, including those actually forming the peak, 
are still linear and Gaussian. 
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FIG. 1. Upper panel: The mixed real-momentum space three- 
point function of equation (4) (solid line) and the perturba¬ 
tion theory result (dot-dashed line) as a function of r. Both 
curves are obtained for q = 0.03 frMpc - , and are normalized 
by Pn n {keq)^{2n/k eq )■ Lower panel: The comparison between 
the two results when the background (calculated from the fea¬ 
tureless power spectrum) is subtracted. 


with S(q, t ) to get 

(4) 

- 2Pn n ((7, t) sin ^ • \7£g(x,t), 

where P\ ln (q,t) is the linear matter power spectrum. 2 

Given that terms of order PimS.g have been neglected 
from the r.h.s., one must ask how accurate is the above 
approximation. Realistically, contains also a smooth 
background. Hence, the approximation is valid so long 
as 


where is the “wiggle” component. In our Universe, 
~ 0-8, the peak location £bao ~ 100/i -1 Mpc, and 
the width a w 10/i -1 Mpc. Therefore, the corrections are 
of order 10-20%, but become larger at the nodes of the 
sine and as q —> 2ircr ~ 1 . 

Note, however, that what is more essential in the above 
derivation is the breakdown of scale invariance character¬ 
ized by a/1 bao “C 1; rather than the actual size of the 
feature. Even for small £™/£ g , the contribution (4), with 
£ g replaced by £)", is distinct though perhaps subleading. 

We tested the above expectations by taking S g to be 
the matter contrast itself, and using the tree-level expres¬ 
sion for the bispectrum in perturbation theory—fig. I. 3 
As seen, subtracting the smooth contribution of the back¬ 
ground results in a much better agreement. 

Squeezed limit bispectrum. — Taking the Fourier 
transform of (4) with respect to x 1 we obtain the 
squeezed limit (q <C k ) momentum space bispectrum: 

( 6(q,t)S g (k-,t)8g(-k + ,t )) 

k (9) 

^ ^2“ P lin {q,t)[Pg(k-,t) - P g (k+,t)], 

where k± = k ± q/2. The above derivation can be gen¬ 
eralized to the case where the fields have different time 
arguments. The result, often called the squeezed limit 
consistency condition (see e.g. [3-5]), reads 


(6(q,t)S g (k-,t 1 )S g (-k + ,t 2 )) ^ 


Piin(q,t) 


r Wi) 

lD(t) 


Pg(k-,t i) 


EM 

D(t ) 


Pg{k+,t 2 ) 


( 10 ) 


where D(t) is the linear growth factor. The k± in the 
arguments of P g are normally approximated by k , which 


3 The tree-level bispectrum is given by 
B(ki, k 2 5 ks) = 2[i^2 {k\ 5 k 2 )P(k±)P(k 2 ) + 2 permutations], (6) 


sin (^y^BAo) y^£ g > (5) 


2 We use finite volume Fourier transformation where the cosine 
mode is related to Fourier modes according to 5 q (t) = [5(q,t) + 
5(—q,t)]/V. However, for convenience, the discrete momentum 
sums are approximated by integrals V f d 3 q/(27r) 3 . Thus, 

(S(q,t)S(q',t)) = P(q)(2ir) 3 6 3 (q + q'), 

with <S 3 (0) = V/(2n) 3 . In what follows, momentum conserva¬ 
tion is always explicitly imposed on momentum-space correla¬ 
tion functions, but the factor (27r) 3 <5 3 (2^ qi) is dropped. Note 
also that the sine modes do not contribute to the relative dis¬ 
placement of pairs located at zby. 


where F n are the usual Standard Perturbation Theory (SPT) 
kernels [1]: 


F 2 (fci, fe 2 ) 


5 

7 


2 (fci -fc 2 ) 2 
7 k^k^ 


+ ^(fei -k 2 ) 




(7) 


For simplicity, the plots are made using the BBKS power spec¬ 
trum [2] modified to account for BAO wiggles 


P(k) = P B BKS«(1 + T W (k/ke q.)) , (8) 

where fc e q = 0.01/iMpc -1 is the equality scale, and the transfer 
function T w (x) is given by 


T w [x) = asm(fx)W(x, #max)(l - W(x,: r m i n )) , 

where W(x,x o) = exp(— x 2 /xq). The parameters are chosen to 
reproduce the observed BAO peak: a = 0.05, / = 1, rr ma x = 30, 
and a: m i n = 3. 
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is valid in the q -A 0 limit: the difference 

P g {k±,t ) - -Pg(M) = ±\<1 ■ V fe P g (M) (11) 

results in an O(q 0 ) term in (10) that is comparable to 
other dynamical effects of the long mode. This has led to 
the conclusion that the 1/q contribution to the squeezed 
limit bispectrum vanishes at equal times. (The 1 /q piece 
survives in unequal time correlations. Measuring unequal 
time correlations is equivalent to watching the galaxies 
as they fall in the long wavelength gravitational field. 
Unfortunately, this is practically impossible.) 

However, the above reasoning does not necessarily hold 
when considering squeezed triangles with small but finite 

q. In the presence of the acoustic feature, P g (k , t) has an 

oscillatory component with period , which can be 

smaller than q. In this regime, the approximation (11) 
is invalid and the difference is proportional to the power 
P™(k,t) in the acoustic peak - the Fourier transform of 
t; g after the subtraction of a smooth background. The 

r. h.s. of (9) now reads 

2P lin (q,t)sm £ B ^g 2 P g ( fc > *)» ( 12 ) 

where To derive this expression, we have 

used the fact that the Fourier transform of a sharp fea¬ 
ture is generically a fast oscillating piece times a smooth 
envelope; equation (8) is an example. For q(BAO 1 the 
result is enhanced by a factor of k/q , 4 

Nevertheless, compared to other dynamical contribu¬ 
tions, expression (12) is suppressed by P™(k)/P g (k). In 
the case of the initial matter power spectrum this ratio 
has support for fc^BAO < 100 and reaches a maximum 
of approximately 0.05 at Mbao ~ 10. The overall result 
turns out to be a subdominant component of the full 
momentum space bispectrum, essentially because most 
of the power at high k comes from short distance cor¬ 
relations £, g (x ~ 2n/k) rather than the acoustic feature. 
A comparison with the tree-level matter bispectrum is 
shown in fig. 2. As seen, once the smooth background 
is subtracted, what remains is well approximated, in the 
squeezed limit, by the universal result (12). 

III. BAO SPREAD AND RECONSTRUCTION 

Intuitively, the above result describes how galaxy pairs, 
which are more likely to be found at distance Ibao, are 


4 In the case of higher point correlation functions (9) generalizes 
to 

( S{q,t)S g (k 1 ,t .) • • ■ S g (k n ,t)) ~ Pn n (q, t) 

Yl ^TT ( a 9( fel U) • ■■Sg(\k + q\,t) • - - <5g(fc„, t)> , 

% 

which again scales as 1/q for q > 27r/:^(o. 



FIG. 2. Upper panel: The bispectrum calculated using equa¬ 
tion (12) (solid line) and the tree-level perturbation theory re¬ 
sult (dot-dashed line) as a function ofk, forq = 0.03 /iMpc -1 . 
Both curves are normalized by P/ in (k eq ). Lower panel: The 
same as above with the smooth background subtracted. 

moved to larger or smaller separations in the presence of a 
mode of wavelength longer than cr. When averaged over 
the long modes, these motions lead to the well-known 
spread of the acoustic peak, as will be discussed in the 
rest of the paper. 

For this purpose, it is necessary to keep higher order 
terms in the expansion (2). At second order in relative 
displacement, now caused by the modes q\ and q- 2 , the 
r.h.s. reads 

2 S qi 6 q2 sin sin (^“) did^ g (x,t). (13) 

As before, this is the leading effect of the long mode if x ss 
^bao j and £ g is the correlation function in the absence of 
the q modes. By correlating (13) with two long modes 
one can obtain the double-squeezed four-point correlation 
function. Alternatively, averaging over the long modes 
with q < A <C 27rcr _1 , gives the first correction to the 
observed two-point correlation around the peak: 

ig(r,t) & tgA r >t) + Zg,s(r,t) + T? A £g !S (r,t), (14) 

where r = |a|, prime denotes d/dr , and terms suppressed 
by ct/^bao are neglected. ^ g ^(x,t) - the direct contri¬ 
bution of the long-modes to the correlation function - 
and £ g ,s(x, t ) - that of the short modes in the absence of 
the long modes- are assumed to be isotropic. Note that 
while s contains the full short scale nonlinearities, only 
the leading effect of the long modes on the short modes 
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has been kept in (14). For each q mode, this scales as 
flin(g , )(^BAo/o ') 2 for q < and P\in(q)/(qcr) 2 for 

q > ^bao- The corrections are suppressed by one or 
more powers of a/t bao and qa , respectively. Hence, due 
to the bulk motions, has a broader peak with E^ given 

by 

S A ~ dgPlin(?)[1 - jo(<*4 ao)+ 2j 2 (^bao )], (15) 

where j n is the n th order spherical Bessel function. 

It is easy to perturbatively confirm the above result 
when £ g is taken to be the dark matter correlation: The 
leading contribution of the long wavelength modes to the 
one-loop power spectrum of the peak reads 5 

pw ( k > A) =1 [ A - d3q (g - fc)2 p, („) 

Pi-loop^ > A ) 2 J ( 27 r)3 qi Fhn[q ’ (17) 

[Pun (I* + Q\) + Pfn(\k - <?|) - 2 PZ(k)} • 


For q <C k the expression in the square brackets simplifies 
t0 -4T > 1 ^(fc)sin 2 (qr • ki B ao/2), giving 

Pi-ioo P (k > A) = S 2 k 2 Pf n (k), (18) 

and taking the Fourier transform with respect to k re¬ 
produces (14). 

Note that for any k, our approximation is valid for all 
q k while the above expressions are based on a rigid 
separation of scales above and below A. Of course, in 
reality P™ (k) has support in a large range of momenta, 
roughly (0.05 — 1) /iMpc . Even if a g-mode falls in this 
range, it is still true that its leading effect on higher k 
modes is the mere bulk motion. Therefore, it contributes 
to the peak power through £ g ,L, and at the same time, 
broadens it by dispersing the shorter modes. A better 
estimate of the width can be obtained by including for 
each k the broadening effect of all smaller q modes, i.e. 
by taking A to increase with k. Below, wc will implement 
this idea by taking A = ek, with e <C 1. 

Taking e = 1/2, the above expression (18) predicts an 
effective broadening of « 5.5/i _1 Mpc, where £:* is 
defined by E^/c* = 1. This turns out to be a sizable 
fraction of the actual width of the observed matter cor¬ 
relation function. We compare the theoretical prediction 


5 The full one-loop power spectrum is given by 

d 3 <7 


I 


oo 3 


[6 F 3 (q, -q, k)P lin (k) + 2F$(q, k - <?)Pi in (|fc - q|)]Hin(«) • 

(16) 

For q <C k it reduces to (17). Incidentally, this coincides with 

1 [ 

2 Ja< 


' q<£.k 




as expected from the remark after (13). 



FIG. 3. The acoustic peak in the matter correlation function 
in linear theory (solid), 1-loop perturbation theory (dashed), 
and simulation. 


with the result of an A'-body simulation 6 in fig. 3. It is 
seen that the perturbative treatment has completely de¬ 
formed the shape of the peak. A more accurate descrip¬ 
tion should, therefore, treat the relative motions non- 
perturbatively. 

Infra-red resummation. — We can obtain a formula 
which is valid to all orders in the relative displacement 
Sq/q, by rewriting (2) as (see e.g. [8]) 


(<«( 2 ■*)«,( 2 


exp 


u ^ ik x 
(2tt) 3 

(■ Sg{k,t)Sg(-k,t )) - 


(19) 


As before, this is only relevant in the presence of a fea¬ 
ture. Taking the expectation value over the realizations 
of the q modes, approximating them, as we did so far, as 
being Gaussian, and using (exp(z^)) = exp(— /2) 
for Gaussian variables, we obtain our final expression 
for the dressed two-point correlation function around 


6 We are measuring power spectra and correlation functions in a 
suite of 16 dark matter only simulations, each of which captures 
the evolution of 1024 3 particles in a box of 1500 3 h _3 Mpc 3 . The 
matter density parameter is f2 m = 0.272, the tilt n s = 0.967 and 
the normalization og = 0.81. The leading cosmic variance has 
been divided out, such that the error bars reflect the sub-leading 
cosmic variance. 
















r « ^BAO 
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/ ^^ m e-^{S a {k,t)S a {rh,t)) € . (20) 

To write the exponent in the above form, we have used 
the fact that V 2 ~ d 2 [and therefore k 2 ~ (x • fc) 2 ] up to 
corrections of order u ji bao- In principle, the exponen¬ 
tial factor should only multiply the peak power P™(k), 
though in practice the smooth background atrw £bao is 
insensitive to the presence of this factor since E <C ^bao ■ 
The subscript e on the momentum space expectation 
value on the r.h.s. indicates that it should be evalu¬ 
ated in the absence of modes with momentum q smaller 
than efc, though it contains all short scale nonlinearities. 
Within a perturbative framework, it is possible to include 
dynamical effects of the long modes, as well as their non- 
Gaussianity by writing more complicated expressions (see 
below). 

To get an idea of how well (20) performs, we set 
S g = 6 and approximate the exclusive expectation value 
in the integral first by the linear matter power spectrum, 
and then by the 1-loop perturbation theory result. The 
first approximation underestimates the broadening by ne¬ 
glecting short scale nonlinearities and therefore predicts 
a slightly sharper peak. 

Let us discuss the 1-loop approximation in more de¬ 
tails to see how (20) can be used to improve perturbative 
results. Two points have to be kept in mind: (i) The 
broadening is only relevant for the acoustic peak, hence 
the exponential broadening in (20) multiplies P™(k). (ii) 
Replacing P™{k) with the 1-loop power spectrum double¬ 
counts the effect of the long modes since the 1-loop re¬ 
sult already contains H 2 k k 2 P™ n (k) [c.f. (18)]. Hence in 
this context the infra-red resummed version of the 1-loop 
power spectrum presented in [7] can be simplified and 
written as: 


P(fc) = PC’(fc) 


pnw 


loop 


(fc) 


( 21 ) 


+e 




(1 + E 2 e]s k 2 )P^(k) + e 


-EAfc" 


1 —loop 


(*)> 


where the first line contains just the smooth part of the 
power spectrum. 7 When considering loop integrals with 
large internal momenta, one should allow for the possi¬ 
bility of higher derivative corrections to the dark matter 
equations of motion in an Effective Field Theory (EFT) 
framework [9]. These corrections compensate for the er¬ 
ror made in treating the short-scale modes as a perfect 
fluid. Therefore, the EFT 1-loop power spectrum differs 
from (16) by one such correction: 

Pi-ioop(fc) = Pia(k) + P 22 (k) - 2P 2 fc 2 P lin (fc), (22) 


where R (also known as speed of sound) is chosen to be 
1.8 hr 2 Mpc 2 in order to obtain 1% agreement with the 


7 In practice, P™ j oop can be obtained by substituting P|j T1 (k) with 
its no-wiggle part in the loop integrals (16) since 1. 



FIG. 4. The ratio of various theoretical approximations to the 
power spectrum to the simulation result. Solid: IR-resummed 
(21), short-dashed: 1-parameter 1-loop EFT (22), dot-dashed: 
0-parameter 1-loop EFT (22) with R = 0, and long-dashed: 
linear. The gray shaded region on the IR-resummed EFT 
curve gives the statistical error. 


simulation results up to fc max = 0.3/tMpc^ 1 (see fig. 4). 
This choice is a rough estimate of P, made in order to 
illustrate how the resummation improves matching the 
BAO oscillations for k > 0.1/iMpc -1 . The exact value of 
R is irrelevant for the shape of the acoustic peak. 

The above resummation formula (21) can be straight¬ 
forwardly extended to any order in perturbation theory 
and to higher order statistics such as the bispectrum or 
trispectrum. Note that in this approximation the lead¬ 
ing dynamical effect of the long modes on short modes is 
also taken into account. The comparison between the IR- 
improved power spectrum (21), and the original 1-loop 
result (22) can be seen in fig. 4. The IR-resummation 
clearly reduces the residual wiggles in the EFT predic¬ 
tion and can thus increase the range over which the the¬ 
ory agrees with simulations, as was pointed out in [7]. 

For the correlation function, the broadened acoustic 
peak resulting from the IR-resummed linear and 1-loop 
power spectra is shown together with the initial peak in 
fig. 5. Although the first approximation does not fully 
capture the smoothing of the peak seen in the data, it 
shows that indeed most of the spread is caused by the 
bulk motions. 

Without resummation the 1-loop EFT (or SPT) power 
spectra result in a spurious double-peaked feature at the 
BAO scale similar to the one shown in fig. 3. This is 
because they only include E 2 fc £"(r) while higher deriva¬ 
tive terms l/n\T, 2 k ^ 2n \r) that partially cancel this fea- 
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ture are not absent. The presence of this feature is the 
cause for the common wisdom that SPT does not work 
for the correlation function. As the good performance of 
the IR-resummed EFT proves, the failure is not related 
to the high-fc behavior of the perturbation theory but 
to the missing non-perturbative treatment of motions. 
One can indeed see that the IR-resummed EFT provides 
a good description of the correlation function down to 
10 h ~ x Mpc separations [7]. 

Another feature of fig. 5 that is worth emphasizing is 
the shift of the peak compared to the linear correlation 
function. This shift is expected to be due to corrections 
to £ g of order E 2 ^/£ba0 5 which are smaller than the 
broadening effects by a factor of ct/^bao [10]- They are 
not entirely fixed by symmetries since the cross corre¬ 
lation between a displacement and other nonuniversal 
effects — e.g. arising from living in an over dense re¬ 
gion — caused by a long wavelength mode contributes at 
the same level. Nevertheless, they can be calculated in 
perturbation theory and are included, to leading order, 
in the 1-loop result, which predicts the position of the 
peak reasonably well. On the other hand, the BAO re¬ 
construction schemes, to be discussed below, reproduce 
the original peak by virtue of undoing the displacements 
caused by the long modes which also eliminates the above 
mentioned cross correlations. 

For comparison, we have also plotted in fig. 5 the 
Zel’dovich correlation function, which is known to give 
a relatively accurate description of the BAO spread. We 
will next argue that the success of the Zel’dovich approx¬ 
imation is because it can be formulated as (20). 

Zel’dovich approximation. — The matter correlation 
function can be related to the correlation function of the 
relative displacement As(z) of two points with initial 
(Lagrangian) separation z: 

l+f(*)= J ^L e ik x j d ^ze- lk z (e~ k As(z) Y (23) 

In the Zel’dovich approximation, As is replaced by its 
linear expression, and the above expectation value is triv¬ 
ially expressed in terms of the variance 

A ij (z) = (As i (z)As J (z)) 

d 3 g ^f--Ffcn(g) sin 2 (—) ■ 

Let us define Zel’dovich power spectrum as the result of 
the inner integral in (23) at k ^ 0: 



FIG. 5. Various theoretical approximations to the acoustic 
peak in the correlation function as well as simulation mea¬ 
surements. Solid: linear, dashed: IR-resummed linear, dot- 
dashed: IR-resummed 1-loop, and dotted: Zel’dovich. 


over short modes q > A, and long modes q < A. So we 
have 


A^(z) = Af(z,A) + A£(z,A). (26) 

A Zel’dovich power spectrum in the absence of the long 
modes P z ,s{k, A), where A -C k, can now be defined by 
replacing A 1 - 7 —> Ag in (25). This is the analog of the 
last factor in (20): it contains the full nonlinear effect of 
the short modes in the Zel’dovich approximation, but no 
long inodes whatsoever. 

Consider now the full P z (k). The integral in (25) is 
dominated by z = 0(l/k), and, if k is in the support of 
P!f (k), by z = ±^bac )k + 0(1 /k). The second contribu¬ 
tion is what we called P™(k). Here, A‘/(z) is first of all 
appreciable, and second, it can be approximated to be a 
constant given by its value at z = £bao k to yield 


P™(k) « e -5 A i(^BAofe,A)fc i fc J p^ g ( fe; A ) 


(27) 


e -^p :<s (k,A). 


P z (k) = J ^ze- ik z e~^ Aii(z)k ' kJ , (25) 

which in the presence of the BAO feature contains an 
oscillating component P™(k). This can be approximated 
by the product of a non-smoothed piece times a broad¬ 
ening factor, as in (20): Define Ag(z, A), and A^(z, A) 
by the same integral as in (24), but taken, respectively, 


The second equality holds up to terms suppressed by 
ct/^bao- Replacing A —> ek results in the desired ana¬ 
log of (20). 

Hence, the Zel’dovich approximation, despite being a 
crude model of short scale dynamics, gives an accurate 
description of BAO broadening by taking into account 
the leading displacement caused by all longer wavelength 
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FIG. 6. The same as fig. 1 with q = 0.1 /iMpc -1 . All curves 
are normalized by Pii n (k e q)£,(2n /fc eq )• Upper panel: Without 
high-pass filter. Lower panel: With high-pass filter. In both 
cases the smooth background is subtracted. 

modes on any given scale k. s 

BAO reconstruction. — This naturally leads us to the 
discussion of BAO reconstruction, and its connection to 
the long-short correlations (4) and (12). The BAO re¬ 
construction is a method to reproduce a sharper acoustic 
peak by undoing the bulk motion induced by the long 
wavelength modes, and hence, it is based on the same 
underlying idea that led to our results [13, 14]. Given 


Two alternative approximations have been proposed in the liter¬ 
ature (e.g. [11, 12]) to model the broadening effects: 


P“(fc)« e - E - fc Pft(fc), 

(28) 

P“(fc) we-^ fca J^ n (fc), 

(29) 


where the velocity dispersion ct 2 is given by the same integral 
as in (15) with A = oo, but without the last square brackets. 
The two expressions happen to give similar results for the mat¬ 
ter correlation function, and to be in good agreement with the 
result of simulations. However, we think the agreement in our 
Universe is accidental. The velocity dispersion is missing the fac¬ 
tor sin 2 (q • x/2) in the relative displacement, which suppresses 
the contribution of the super-long modes. Had there been more 
power at large scales, or if fceq^BAO <C 1, (29) and (28) would 
have differed significantly. On the other hand, equation (28) ap¬ 
proximates the short-long effects by the same expression as that 
of the long-short effects. This is not justified by any symmetry 
argument, and is an overestimation in the real universe. (28) 
would predict too much spreading if there was more power in 
small scales. 


that the leading effect of the long mode is a uniform 
acceleration, this procedure roles back part of the time- 
evolution, which as we saw leads to the broadening of the 
BAO peak. 

Operationally, the reconstruction method consists of 
three steps: (i) Choosing a rigid separation A between 
long and short modes, and solving for the linear dis¬ 
placement field produced by the long modes, (ii) Mov¬ 
ing back all points according to this linear displacement 
field (as one would do in the Zel’dovich approximation), 
(iii) Adding back the original smooth field that is largely 
erased by step (ii). In this procedure, the only effect of 
the long modes on the short modes that has been reliably 
taken into account is the linear displacement. Hence, the 
effectiveness of the method seems to be a strong indica¬ 
tion of the validity of (4). But, there are two caveats. 
First, the reconstruction method does not significantly 
affect the smooth part of the correlation function, hence 
it only verifies (4) after background subtraction. 

Second, the threshold A is practically chosen within 
the support of P ™, where as mentioned above, the modes 
both contribute to the peak, and cause it to spread. It is 
natural to suspect equation (4) to become a poor approx¬ 
imation for these q modes, due to their dynamical self¬ 
coupling. On the other hand, the reconstruction method 
would primarily deal with the effect of the q modes on 
higher k modes. Therefore, the effectiveness of recon¬ 
struction implies that even for these relatively larger val¬ 
ues of g, once the contribution of modes below the thresh¬ 
old A is removed from S g , equation (4) should be a good 
approximation. 

To test this expectation in perturbation theory, we in¬ 
sert a high-pass filtered power spectrum P g (k) = (1 — 
W(k,q))P]in(k) into the tree-level matter bispectrum, 
while keeping Pii n (g) unfiltered. The inverse Fourier 
transform with respect to k is then to be compared to 
the r.h.s. of (4), with £ g obtained from the same high- 
pass filtered P g (k). The results are shown in fig. 6, and 
seem to be in moderate agreement. The high-pass fil¬ 
ter effectively picks small laboratories, free-falling in the 
background of the long wavelength mode. 


IV. CONCLUSIONS 

We used the leading Newtonian effect of a long wave¬ 
length matter perturbation Sr to derive approximate for¬ 
mulas for its correlation with the distribution of pairs (4), 
as well as the squeezed limit bispectrum (12), in the pres¬ 
ence of the BAO feature. The derivation is based only on 
two underlying assumptions: first, the equivalence prin¬ 
ciple, by which we imply that no additional (fifth) force is 
universally sourced by material objects, and second, lo¬ 
cal formation of tracers which forbids nontrivial bias with 
respect to locally unobservable quantities such as veloc¬ 
ity and gravitational potential. This requires absence, 
or rather smallness, of primordial local non-Gaussianity. 
Therefore, the result holds beyond the standard pertur- 










bation theory, and apply equally well to biased tracers. In 
the real Universe, it gives the dominant component of the 
real space correlation at x ~ ^baO; but a subdominant - 
though still unique and distinguishable-piece in momen¬ 
tum space. 

Next, we explored the connection with the broaden¬ 
ing of the acoustic peak, where the same universal effect 
but averaged over the long modes is known to account 
for most of the spread in the observed Universe. We de¬ 
rived a formula for the observed correlation function (20), 
which resums the induced motion by the long modes to all 
orders. A simpler way to implement this IR-resummation 
in perturbation theory was proposed, and the result was 
shown to be in good agreement with the numerical re¬ 
sults from a TV-body simulation, and with the Zel’dovich 


approximation. It was shown that the Zel’dovich approx¬ 
imation to the correlation function can be recast into 
the form of our IR-resummed formula (20), which we 
take as the explanation for its success in predicting the 
BAO spread. Finally, we discussed BAO reconstruction 
method as a practical application of the same underlying 
idea. 
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